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■ Abstract. We present an example of quantum process tomography (QPT) performed 
\^ ! on a single solid state qubit. The qubit used is two energy levels of the triplet state 

in the Nitrogen- Vacancy defect in Diamond. Quantum process tomography is applied 
to a qubit which has been allowed to decohere for three different time periods. In 
each case the process is found in terms of the x matrix representation and the affine 

■ map representation. The discrepancy between experimentally estimated process and 
the closest physically valid process is noted. The results of QPT performed after three 
different decoherence times are used to find the error generators, or Lindblad operators, 
for the system, using the technique introduced by Boulant et al 6 . 
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1. Introduction 

Quantum process tomography (QPT) is a method of experimentally determining the 
unknown dynamics of a quantum system. This is crucial as a method of checking the 
functionality of would-be quantum information processing devices. Using knowledge 
obtained from QPT applied to such a device, one can ideally locate and possibly rectify 
any sources of errors or decoherence. The standard QPT technique, which involves the 
use of multiple test states, was first developed by Nielsen and Chuang jT] and Poyatos, 
Cirac and ZoUer pj. A different technique, exploiting entangled states, was subsequently 
proposed by Leung [3j and D'Ariano and Lo Presti [4J. More recently, methods 
have been developed to determine a process from a tomographically incomplete set of 
measurements [5J , as well as techniques to ascertain the master equation describing time 
evolution of the system jEllZI- QPT has been performed in liquid NMR implementations 
|Hl El QUI , numerous optical systems El El El E| , atoms in an optical lattice 
fT7j and in bulk solid state NMR JH|- To the best of our knowledge, this work contains 
the first example of QPT performed on a single solid state qubit. 
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2. The Nitrogen Vacancy Defect in Diamond 

The nitrogen vacancy defect (or N-V center) is a naturally occurring defect in diamond 
with nitrogen impurities. It can be manufactured in type IB synthetic diamond by 
irradiation and subsequent annealing at temperatures above 550 C; radiation damage 
causes vacancies in the diamond lattice and annealing leads to migration of vacancies 
towards the nitrogen atoms (see figure 1(a)). The N-V center can also be produced in 
type IIA diamonds by ion implantation |TI?1 OHj . 




Figure 1. (a) Structure of the N-V defect (N=Nitrogen, V= Vacancy and C=Carbon) 
(b) Energy level scheme of N-V defect in a weak magnetic field. and are excited 
and ground triplet states respectively. The = and = 1 states are used as a 
logical qubit 



The energy level scheme of the N-V center is shown in figure Q^b) . It consists of a 
triplet ground state and a triplet excited state with a metastable singlet state ^A. 
At zero magnetic field the ground triplet state energy levels are split into degenerate 
777.5 = ±1 sublevels and the = sublevel. This is due to the magnetic dipole-dipole 
interaction between the 2 unpaired electrons. This dipole-dipole interaction can be 
described by a Hamiltonian term = S • D • S where D is the zero-field splitting (or 
fine structure) tensor and S = {Sx^ Sy^ S^). Using standard techniques ^22j this dipolar 
term can be rewritten as 

^ +E{Sl + S'^), (1) 



Hn = D 



St - -S{S + 1) 



where D and E are zero-field splitting parameters which are dependent on the symmetry 
of the molecule. The axial symmetry of the N-V center means that = MHz, while 
D for this system has previously been characterized jSHll^ as D = 2880 MHz. The end 
result is that, at zero magnetic field, the = ±1 sublevels are degenerate and at an 
energy 2880 MHz greater than that of the = sublevel. With an applied magnetic 
field aligned along the z axis of the molecular system (i.e. B = (S^,0,0)) the total 
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(2) 



where the gyromagnetic ratio, gf3^ is 2.8025 MHz/Gauss. This Hamiltonian leads to a 
|0) ^ |1) transition of energy 

A^; = (2, 880 + 2.80255^)MHz (3) 

where is the applied magnetic field (figure [21(a)), measured in Gauss. 
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Figure 2. (a) Energy level schematic of spin S=l triplet in an applied magnetic 
field. Because of the axial symmetry of the N-V defect, the ms = ±1 sublevels remain 
degenerate in zero field and are split in non-zero magnetic field, (b) State readout 
is via Rabi oscillations of fluorescence intensity. The circles represent experimentally 
measured data points. 



The fiuorescence intensity corresponding to the = transition between the 
ground {^A) and excited {^E) triplet states is strongest. Optically detected magnetic 
resonance (ODMR) is used to perform state readout [2Sl EE! EZI , i-e. the population 
of the rris = sublevel is identified by the amount of measured fiuorescence (figure 
[21(b)). Manipulation of the qubit state is via standard electron spin resonance (ESR) 
techniques, using microwave pulses resonant with the |0) ^ |1) transition 

The spin longitudinal relaxation time, Ti, is on the order of milliseconds at room 
temperature [2H| (relaxation time of the order of seconds is expected at low temperature) . 
Transverse relaxation times (or decoherence times), T2, of up to 60 microseconds have 
been reported [29j, for samples with low nitrogen concentration. 

In the N-V center, initialization into the = state is achieved by optical 
pumping. Optically induced spin polarization is thought to be related to spin-selective 
intersystem crossing from the photoexcited triplet state to the metastable singlet 
state (Figure njb)) i.e there is an irreversible transition between ^E ^A. The spin 
polarization achieved by this pumping corresponds to at least 70% population in the 
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rris = sublevel at room temperature. Instead of a pure state, corresponding to 100% 
polarization, we have what is known as a pseudopure state which takes the general form 

Ppseudo= + (4) 

where n is the number of qubits and < a < 1 is related to the amount of polarization. 

In the case of the N-V center, the 70% population of the = sublevel, after 
optical pumping, corresponds to TT{ppseudo\0){0\) = .7 and so 

Ppseudo=-Q + .4|^)(V^|. (5) 

It is hoped that, in future experiments, the use of projective readouts, which are 
possible at low temperatures, will lead to increased polarizations. 

A controlled two-qubit quantum gate in which the vacancy electron spin is hyperfine 
coupled to a nearby C^^ nucleus, has recently been performed using this system jHOj- 



3. Standard QPT 

3.1. QPT Technique 

We briefly review the standard QPT technique: 

• Prepare a complete basis of input states pi . • • Pd? (e.g pj = |'0j)('0jU'0j) = 

{|0), |1), ^(|0) + |1)), ^(|0) + i|l))}, for a single qubit) 

• Apply the unknown process S to each member pj 

• Reconstruct the output states £{pj) using quantum state tomography 

While each of the reconstructed £{pj) will be physically valid density matrices (by 
appropriate use of quantum state tomography) it is possible that the process £ describes 
is unphysical e.g. the £{pj) are mutually inconsistent or £ is not completely positive. 
Experimental noise and a finite number of measurements to determine 
expectation values can result in a reconstructed £ which is not physically 
valid. 

If it is possible to write a process in the Kraus representation then we can be sure 
it represents a physically valid process (one that is described by a completely positive 
map). We can specify an unknown process £ by experimentally determining the Kraus 
operators [Si] {Ei} which describe it 

d^ 

£{p) = Y.^^P^l' (6) 

or, using a fixed (complete) basis of operators {Ai\ 

d^ 

^(P)= XmnAmpAl. (7) 

m,n=l 
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Xmn is a matrix of coefficients which completely describes the process S and is 
positive Hermitean by construction. The trace preserving constraint J^i^l^i = ^ 
becomes J2m,n XmnAlAm = I. 

• From the measurement results {^(pi) . . .£{pd'^)}^ Xjk can be determined, given the 

relation S{pj) = XjkPk- 

• In order to determine x from the matrix A one operates on A with the pseudoinverse 

of /?, where /? is derived theoretically from the relation A^pjAl = Pjl^Pk- 

• Using this last relation and ([Zj) we can see Xij = "^j^n f^ij^Xmn and so inverting /? 

gives us X 5 as required. 

To derive {Ei} from x we first diagonalize x with a unitary W 

Xmn = ^ Umxdk^klU^l (8) 
k,l 

where di are the eigenvalues of x- The Kraus operators can then be obtained by 

j 

Note that this procedure only works if rf^ > which follows from the positivity 
(semidefinite) of x- It is this property which we use to check the physicality of a process 
£] if the X matrix reconstructed from experimental data has negative eigenvalue(s) 
then this indicates that noise and/or finitely sampled expectation values has caused 
the output data to infer an unphysical process. To overcome this problem, a physical 
matrix x is found which is as close as possible to the original x some sense. Specifically, 
we used a technique analogous to MLE for state estimation j22 E| by minimising a 
deviation function, A(t), incorporating a general parameterization for a positive x • 



= Y \Xmn{t) - Xmnf + X 
m,n=l 



J2 Xmn{t)AlA^-I 



=1 



(10) 



where A is a Lagrange multiplier (which ensures a trace-preserving process) and 

x = THt)m (11) 

where T{t) is a rf^ x rf^ complex, lower triangular matrix with d^ real parameters, t{i). 
In the case of a single qubit we have: 

m 



( ^(1) 








\ 


t(5) + it(6) 


t{2) 








t(ll) +it(12) 


t{7) + it{8) 


m 





\ t(15) +it(16) 


t{13)+it{U) 


t{9) + it{10) 


^(4) / 



(12) 



The algorithm used to minimise A(t) in ffTUj) was the Nelder-Mead simplex algorithm 
as implemented in Matlab @. 
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3.2. Avoiding Local Minima 

As is often the case with numerical minimization problems, (pT)|) typically contains 
numerous local minima. Our preliminary investigations suggested that, given a random 
initial point in parameter space, the algorithm would often fail to find the global 
minimum. As such, a good starting point, in f[T^ was deemed necessary if 

the results were to be meaningful. In order to obtain we used a technique based 

on principal component analysis [SSI- If experimentally determined matrix x is not 
positive semidefinite, one can "filter" it by setting any (presumably quite small) negative 
eigenvalues to zero. Specifically, if we decompose x 

X = UDU\ (13) 

where D is a diagonal matrix containing the eigenvalues of x ? then we can construct a 
similar but positive matrix via 

X* = UD^U\ (14) 

where is identical to D except that any negative eigenvalues have been set equal to 
zero. From x^ we can extract good initial parameters t^{i) by performing a Cholesky 
decomposition (e.g. using the in-built "Choi" function in Matlab ®). 



3.3. Process Visualisation 



Any one qubit state can be parameterized by a so-called Bloch vector r. Equivalently, 
one can explicitly include the unit coefficient of the identity basis component and 
parameterise the qubit state by a 4-vector (1, r^^, r^, r^)^: 

/ 1 \ 



p 



-(I + r-a) 



(15) 



where + + < 1 [Sj. In this basis any linear trace-preserving evolution 8 takes 
the aflSne form 

/ 1 \ 






FEE 

^xx ^yx 

-^Xy Eqiqi E 



yy 



'zy 



f E E E 

^7. ^T,7. yz ^ '■ 



(16) 



and the process applied to an arbitrary input state is an affine map on r, (although it 
remains a linear map on p) : 



(17) 






or, more explicitly. 



rf = Er + t. 



(18) 
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Identifying the surface of the unit Bloch sphere as the set of all possible 
input states, we can visualize a process S by its action on the sphere. The 

3x3 matrix E is responsible for deformation and rotation of the Bloch sphere, while t 
denotes displacement from r = (0,0,0). 

We can gain crude but immediate insight into the nature of a single qubit process 
by using this visualization technique. Bloch vectors with |r| > 1 do not correspond to 
any valid density matrix. Hence, when a process is depicted this way, any protrusions 
of the output ellipsoid outside the unit Bloch sphere mean that the process is not trace 
preserving. 

Complete posit ivity also places constraints on physically allowable output ellipsoids. 
For a more detailed analysis of this problem we refer to [SSI EH E^- The crucial point to 
note is that seemingly innocuous ellipsoids (i.e. contained within the unit Bloch sphere 
and not overly deformed) can represent a process which is not physically valid. The 
variety of allowable ellipsoids is particularly reduced when the translation vector t in 
is non-zero. 

4. Experimental Method 

The QPT experiment was performed at room temperature using diamond nanocrystals 
obtained from type lb synthetic diamond. Diamond nanocrystals were spin coated 
on a glass substrate, and single nanocrystals were observed with a homebuilt sample- 
scanning confocal microscope. In order to ensure the presence of a single defect in the 
laser focus, the second-order coherence was measured using Hanbury-Brown and Twiss 
interferometer and the contrast of the antibunching depth was determined. In order to 
perform many repetitions (to obtain good expectation values) in a reasonable amount 
of time, a sample with a relatively short coherence time (2/xs) was chosen. Microwaves 
were coupled to the sample by a ESR microresonator connected to a 40W travelling 
wave tube amplifier. A magnetic field {B = (0, 0, 5^), = 200 Gauss) was applied to 
the defect in order to remove degeneracy of the | ± 1) energy levels. 

Rabi oscillations of fiuoresence intensity were used to obtain expectation values. A 
reference nutation was initially taken in order to normalize the measurement nutations 
(i.e. set 0,1 levels for expectation values for the pseudopure initial state). In addition 
this reference oscillation was used to derive the microwave pulse time required to perform 
I and TT rotations of the Bloch vector. A complete basis of four input states was then 
prepared using microwave pulses resonant with qubit transitions. The po state was 
obtained directly by optical pumping and three remaining input states pi, p2 and ps 
were obtained by application of suitable tt or | pulses. Each of these states was left to 
decohere for a series of time intervals r. As a last step, measurements of the diagonal 
and off-diagonal elements of the density matrix were performed. 

Estimates of the density matrix elements were extracted from experimental data 
(Rabi oscillations) using the maximum entropy (MaxEnt) technique jHSI- This method 
returns a physically valid density matrix which satisfies, as closely as possible, the 
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expectation values of measured observables. In cases where only incomplete knowledge 
of the output state is known, an additional constraint is used; the reconstructed state 
must also have the maximum allowable von Neumann entropy. 

5. Quantum process tomography experiment 

In summary the QPT experiment was performed as follows: 
• The N-V center was optically pumped to the initial pseudopure state 



Hereafter we neglect the identity component and treat this initial state as p = |0)(0|. 

• Following this, each of the following input states pi . . . (pj = = 

{|0), -^(lO) + |1)), -^(lO) + was created by application of suitable tt or | 

pulses. 

• Each input state was then left to decohere, solely by interaction with its surroundings 

(i.e. no microwave (ESR) or optical (ODMR) radiation was applied to the system), 
for a decoherence time, r. 

• After T had elapsed, expectation values for dx^cfy and dz were obtained by ODMR. 

• Using < dx >^< (fy > and < dz > the final (output) states £{pi)^ after the 

decoherence, were reconstructed via the MaxEnt technique. 

• Using these output states, the process for each time period was reconstructed using 

the technique described above i.e. minimization of ([Tn|) 

The processes, both from raw data and reconstructed (i.e. completley positive), 
are depicted in figures El to El 

5.1. Quantification of Unphysicality 

The Jamiolkowski isomorphism [32j maps a quantum operation ^, acting on Hd^ to a 
quantum state pg in Hd'^ via 



where {\j)} is some orthonormal basis set. 

When the process £ is physical, one then obtains a physically valid ps which can be 
compared to the ideal process (converted to pid) using distance measures on quantum 
states The trace distance between density matrices pid and ps is 




(19) 




(20) 



(21) 



D{pid,pe) = -tr I Pid - p£ 



(22) 



Quantum process tomography of a single solid state qubit 



9 




Figure 3. 20ns results: (a) Process from experimental data, x- (b) Process obtained 
from physically valid x- 





Process Matrix £j[ 


Experimental 




( 1.0000 \ 
0.0626 0.6552 -0.0225 -0.1198 
0.0448 0.0287 0.7309 -0.0226 

^ 0.0138 -0.0143 0.0878 0.9843 ) 




Reconstructed 




( 1.0000 ^ \ 
0.0532 0.6798 -0.0312 -0.1093 
0.0420 0.0206 0.7051 -0.0227 

^ 0.0070 0.0001 0.0916 0.9410 y 





where 



X \= Vxtx. 



(23) 
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Figure 4. 40ns results: (a) Process from experimental data, x- (b) Process obtained 
from physically valid x- 





Process Matrix £j[ 


Experimental 




/ 1.0000 \ 
-0.0344 0.6378 0.1660 -0.0261 
0.0683 -0.0728 0.6872 0.0499 
I 0.0277 0.0536 0.0827 0.9614 J 




Reconstructed 




/ 1.0000 1 
-0.0307 0.6477 0.1443 -0.0266 
0.0638 -0.0944 0.6773 0.0407 
^ -0.0004 0.0532 0.0824 0.9320 y 





Similarly one can define the Fidelity: 

F{pid,p£) = tr 



PsVPid J ■ 



(24) 
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Figure 5. 80ns results: (a) Process from experimental data, x- (b) Process obtained 
from physically valid x- 





Process Matrix £j[ 


Experimental 




/ 1.0000 \ 
0.0442 0.2359 -0.1001 -0.0757 
0.0791 -0.0947 0.3770 -0.1163 

I 0.0434 0.0487 -0.0461 0.9554 J 




Reconstructed 




/ 1.0000 1 
0.0503 0.2726 -0.0290 -0.0752 
0.0834 -0.0357 0.3291 -0.1038 

^ 0.0302 0.0457 -0.0495 0.8984 y 





Using this definition we define the Bures Metric 

B{pid, Ps) = \/2-2Vi^(Ad,P^ (25) 
and the C Metric 

C{pid, Ps) = ^l-F{p,a,Ps). (26) 
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An unphysical process, however, can lead to an unphysical ps^ possibly resulting in a 
process fidelity which is greater than one. The application of the preceding fidelity- 
based distance measures can, therefore, produce nonsensical results. In such cases it is 
necessary to use other techniques in order to estimate the disparity between x X 
for example. If one defines X = x ~ then possible measures are the matrix p-norms 
{p = 1, 2, oc) of X and the Frobenius norm of X {\\X\\Fro) as well as the trace distance 

As stated above, these quantities gives a measure of how well x describes the 
experimental results. The difference matrix X varies depending on the basis of operators 
chosen for x- We propose that if a standard basis is chosen then the norms of X, 
defined above, provide a method of "benchmarking" the quality of all such quantum 
process tomography experiments, regardless of implementation. Table Q] describes the 
discrepancies between experimental and reconstructed processes. Here we have chosen 
the normal basis {Ak}{Adi+j^i = \i){j\) with which to construct x x • This 
particular choice of basis means that x is equal to the Choi matrix C jMl lE^ 

d-i 

X = C=J2^(\')0\)^\^){jl (27) 

and proportional to the characteristic state ps for the process: 

X = dps = d{£®Id)m{^\. (28) 



Decoherence 
time 


P^llp=i 






Dpj-o 


20ns 


0.0660 


0.0525 


0.0636 


0.0262 


40ns 


0.0581 


0.0427 


0.0529 


0.0213 


80ns 


0.1141 


0.1075 


0.1276 


0.0538 



Table 1. The disparity between the estimated process from the experimental 
measurements (x) and the "nearest" physicahy valid process (x). X is defined as 
X = X-X' 



6. Markovian Process Tomography 

Standard QPT utilises a "black box" approach to studying dynamics. It predicts 
the resultant output states given arbitrary initial states, but fails to describe the 
path taken through state space, over time, from initial to final state. If we are 
prepared to make certain assumptions about the relationship between system and 
environment, the Born and Markov approximations, then we can construct a Markovian 
quantum master equation which describes the time evolution of the system over the 
time studied. In the area of quantum information processing, understanding open 
system dynamics is important in studying quantum noise processes, designing error 
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correcting codes ^01 and locating decoherence free subspaces [¥T]. We discuss a method 
of experimentally reconstructing a master equation by using estimates of the state of 
the system at a number of different timepoints. This particular technique was developed 
and implemented by Boulant, Havel, Pravia and Cory [H], in the context of liquid NMR 
quantum computation. The method is different to standard QPT techniques insofar as 
it assumes Markovian dynamics and, consequentially, the reconstructed equation has 
separate terms which describe the unitary and non-unitary aspects of the evolution. We 
then proceed to apply this technique to the N-V center, in an effort to better identify 
the decoherence processes it suffers. This is the first time such an analysis has been 
performed on a single solid state qubit. 

By invoking the Born and Markov approximations one may derive the so-called 
Lindblad Master Equation commonly expressed in a number of different but 
equivalent ways, including: 



dp 
di 

dp 
di 

dp 
di 



i[H,p] + 



-Y 

k=l 



(2LkpLl - LlLkP - pLlLk^ , 



k=i 



(29) 



(30) 



(31) 



The Lindblad master equation dHOI) is the diagonal form of the GKS master equation 
derived by Gorini, Kossakowski and Sudarshan 



dp 
di 



(32) 



where aa(3 is a {dP — 1) x [d^ — 1) Hermitean matrix and {F^,} is a linear basis of traceless 
operators on density matrices. If the matrix of coefficients a^,^ (sometimes called the 
GKS matrix ^Hl) is positive, then the process described by f|H^ is completely positive. 
Diagonalising a^p in f|H^ leads to the Linblad form f|Hn|) . 

We can clearly see one advantage gained by making the Born-Markov 
approximations - the resulting master equation can be separated into components which 
cause unitary and non-unitary evolution of the system: 



Unitary Evolution : 

Non — Unitary Evolution : 



t[H,p], 



-y 

2 ^ 

k=l 



LkPi Lf, 



+ 



Lk,pLl 



(33) 
(34) 



The operators are Lindblad operators and they describe the decoherence processes. 

We now describe a technique for estimating the generator >C of a Markovian process. 
We transform to a Liouville space basis, as in [6j, where density matrices become 
column vectors (denoted and dynamical maps become d^xd^ matrix superoperators 
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(denoted by " ) : 

d\p) 



jC\p). (35) 



dt ^' 

The matrix superpropagator V{t) is then defined by exponentiation: 

ipjt)) = e %m = ht)\pm- (36) 

Henceforth we adopt a similar notation to [01 : 



iH,P] + ^J2{ h/^' 4] + [Lk, pLl] ) (37) 
k=i 



m\pi - n\pi (38) 
in + n\\p) (39) 



dt 

becomes 
dt 



We call H the Hamiltonian superoperator and TZ the relaxation superoperator. 
Exponentiating the generator gives us the superpropagator 

V{t) = e-(^^+^)^ (40) 
We can manipulate (^01) to isolate the relaxation superoperator: 

7^ = -in - hn (V^ (41) 

An alternative method of deriving TZ from estimates of V and TC is to estimate the 
derivative at t = of 



^yn^-[in+iz)t^,^.n_ (42) 

Here we need to invoke the symmetric Baker-Campbell-Hausdorff formula, for arbitrary 
operators A and 5, 

Making the identification 

A = iH, (44) 
B= -[iH + n), (45) 

we get 

ei^^e-^^'-^yei'^ = e"*^ + 0{f). (46) 

Clearly, differentiating f|inj) at t = gives us 7^, as required. In order to obtain 
this diflFerential we used a numerical differentiation technique - Richardson extrapolation 
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To use this technique we require estimates of the propagator, = V{tm)^ at time 
points tjn = 2^ti. This alone, however, is not a robust method for estimating TZ. It 
magnifies the noise present in the estimate of Vi [6j. More importantly, there are no 
constraints on the physicality (i.e. complete positivity) of the generator. This estimate 
of TZ which we will call TZre is instead used as the starting point for a constrained fit 
to the experimental data. 

In order to search for a physical process which was close to the measured data we 
used a parameterisation based on the Gorini Kossakowski Sudarshan master equation: 



C = -i[H,p] + - aaf3 



(47) 



Gorini et al. have proven ^42j that C is the generator of a Markovian semigroup on 
d-dimensional Hilbert space if and only if it can be written in the form (jTTj) . where 
the GKS matrix aai3 is a {(P — 1) x (cP — 1) positive (semidefinite) matrix and {Fa} 
is a hnear basis of traceless operators on p. We will exploit the requirement that Gap 
be positive in a way that is completely analogous to a technique we used in standard 
QPT. During process reconstruction we enforced positivity of x, and, by extension, the 
complete positivity of S{p). Similarly we will enforce complete positivity of the 

decoherence process TZ by constraining the GKS matrix a^/? to be positive. 

For one qubit we have a parameterisation in terms of 9 real numbers x(i): 

d?-i 



FaP,F, 



+ 



Fa,pF, 



(48) 



a,l3=l 



where 





x{2) 








\ 



(49) 



/ x{l) 

x{4) + ix{5) 
^x{8) + ix{9) x{6)+ix{7) x{3) J 

There is considerable freedom in which basis {Fa} to use but we chose to use 
generators for SU(d), rescaled to be trace orthonormal: 

J_( 1 
V2 \ 1 



Fi 



Fo = 



F.= 




(50) 



^2 I -1 



(51) 



(52) 



A numerical search, using Matlab's® fminsearch algorithm, was used to find the 
minimum of 



M 



A{n{x)) = E 



m=l 



H+n{x)\tm ^ 

6 ^ ^ / r. 



(53) 
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where Vm are the experimentaUy determined propagators, for evolution time i.e. 



V = 









( \ 






( \ 






\pl{tm)) 




\p2{tm)) 


|p3(tm)) 
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where 



|pi(0)) = 



/ 1 \ 






|P2(0)) = 



1 





|P3(0)) = 





1 



|P4(0)) = 



(54) 








We label the GKS matrix which most closely fits the data, while remaining physical, 
as (see figures IHl and [Zj) . To derive the Lindblad operators {Lj.} from we first 
diagonalise with a unitary [/''', 

= ^ UmxdxSxyU^y, (56) 



xy 



where di are the eigenvalues of aafs- The Lindblad operators can then be obtained by 
L, = ^,J2Uj,F,. (57) 

These are depicted in figures IHl to fTUl The expectation values < dx >,< (Jy > and < dz > 
for both the experimental and the reconstructed (Markovian, completely positive) 
process are compared in figure El It is natural now to ask given the estimations for L^, 
whether one can deduce any information regarding the physical decoherence processes 
responsible. Unfortunately, this is not generally possible. The most general Markovian 
master equation for a qubit can be fully described by three Linblad operators. However 
this compressed description often will not correspond to the most appropriate physical 
description of the decoherence. It can be argued however, that in most situations, 
e.g. decoherence free subspace investigation, decoupling pulse generation etc., all the 
relevant information is contained in the compressed Lindblad description. 

The relative contribution of each Lindblad operator, L^, to the overall decoherence 
process can be quantified using their Frobenius norms: 



Relative Contribution(Li) 



|Li" 



Fro 



-■^jlFro 



(58) 



In summary the technique we used for ascertaining the relaxation superoperator TZ 
is as follows: 

• Initialisation via optical pumping and appropriate ESR pulses to create a complete 

basis of input states pi. 

• Each Pi was left to decohere for a sequence of times = 2^ti i.e. {20ns, 40ns, 

80ns}. 
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• State tomography was used to determine the output state Pi{tjn) after each 

decoherence time t^- The results were used to construct the matrix 
superpropagators Vm- 

• Using Richardson extrapolation, an estimate of the generator, 1Zre-> was extracted 

from the measured Vm and used as a starting point for the fitting procedure (|SS|)- 

• The minimisation of (|53|) produced a Markovian generator, TZ^ which was physically 

valid and best fit the measured data at a sequence of timepoints t^- 

• The GKS matrix which minimised (|SS|), gl^(3) was then diagonalised in order to find 

the Lindblad operators {Lk}- 

6.1. Markovian process tomography results 




Figure 6. GKS matrix ao^f^: the initial starting point for minimization of (a) 
Real components of aoc(3 (b) Imaginary components of aoc(3- 




Figure 7. GKS matrix a^^'- the matrix which minimizes (|53j) (a) Real components of 
(b) Imaginary components of o^. 
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Figure 8. Lindblad operator Li, Relative contribution ^ 94.4%: (a) Real components 
of Li. (b) Imaginary components of Li. 



0.0829 - O.OOOOi -0.0056 - 0.007H 
-0.0011 + 0.010H -0.0829 + O.OOOOi 




7. Summary 

We have presented the first quantum process tomographic analysis of an individual single 
solid-state qubit, a Nitrogen Vacancy centre in Diamond. This analysis is only possible 
due to the enormous advances made in recent years in single-molecule spectroscopy 
[2Sl EE! EZI , where the resultant ODMR technique provides us here with high-fidelity 
single-qubit readout. As experimental refinements and technological advances improve 
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the readout and control of this system, QPT wiU become an even more important 
diagnostic tool. For example, improvement of the = polarisation by optical 
pumping and increased accuracy of rotations by the use of composite pulse sequences, are 
both feasible short-term goals. If states can be prepared, controlled and read out with 
very high accuracy, then any deviations from Markovian dynamics (e.g. the disparity 
in expectation values, between experimental and reconstructed processes apparent in 
figure fTT|) cannot be dismissed as noise and the nature of the non-Markovian environment 
should be investigated. If two qubit gates can be performed in this system then QPT can 
be used to verify the robustness of encoded information in decoherence free subspaces. 
As quantum devices develop and increase in size, the task of "debugging" the device, or 
actively identifying the noise present in the device, will pose significant challenges. The 
work presented here represents an initial step towards the testing of quantum devices 
in solid-state. 
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Figure 11. Comparison of experimental expectation values at timepoints tm with 
expectation values predicted by Markovian master equation. The time is plotted on a 
logarithmic scale for clarity. 
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